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ABSTRACT 

We investigate the stellar composition of bulges and elliptical galaxies as predicted by 
the CDM paradigm using semi-analytical modelling. We argue that spheroid stars are 
built up of two main components, merger and quiescent, according to the origin of the 
stars. The merger component is formed during major mergers by gas driven to the 
centre, while the quiescent component is formed in gaseous discs and added later to 
the spheroid during major mergers. Galaxies more massive than Mc = 3 x 10 10 M Q 
have on average only a 15% merger component in their spheroids, while smaller galax- 
ies can have up to 30%. The merger component increases with redshift due to mergers 
involving more gas. However we do not find mergers with gas fraction above ~ 40% of 
the remnants mass. Generally the gas fraction is a decreasing function of the redshift 
at which the merger occurs and the mass of the remnant, with more massive remnants 
having smaller gas fraction and hence smaller merger components. This trend is in- 
dependent of the environment of the galaxy with the only impact of the environment 
being that galaxies less massive than Mc have slightly larger merger components in 
dense environments. The fraction of stars in bulges for galaxies more massive than 
Mc is larger than 50%. We find that the majority of stars in galaxies more massive 
than Mc reside within bulges and ellipticals independent of redshift and that the frac- 
tion increases with redshift. The most massive galaxies at each redshift are elliptical 
galaxies. 
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1 INTRODUCTION 

Bright elliptical galaxies are old stellar systems which 
formed the bulk of th eir stars at high re d shift within a 
short time scale (e.g. ICimatti et al.l l2004t iThomas et alJ 
2005). However, they are not a coeval family of objects 
as they for example divide into luminous slow rotating 
pressure supported ellipticals, c alled boxy, and fas t rotat- 
ing less luminous ellipticals fe.g. lBender et al"lll993l) . called 
disky. Another devision of elliptical galaxies suggests that 
they can be divided into core and power-law galaxies ac- 
cording to the steepness of t heir inner surface brightness 
profile fe.g. lLauer et alJ IT995'l. Furthermore, 2-dimensional 
integral- field spectographs like SAURON or OASIS, reveal 
kine matically decoupled co res in many elliptical galaxies 
(e.g. iMcDermid et alJl2005t) . 

There are two competing scenarios for the formation 
of elliptical galaxies. One assumes that they are formed at 
high redshift during a highly efficient starburs t followed by 
passive evolution of the stellar component fe.g. lLarsonll974 : 
IChiosi fc Carraroll2002ft . and the other assumes the forma- 
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tion of an elliptical galaxy in the merger of two galaxies 
(e.g. iToomre fc Toomrelll972T) . In this paper, we will test 
the predictions made by the latter model. 

Since the early proposal of IToomre fc Toomrd <ll972T) 
that early type galaxies form as the result of a binary merger 
between two spiral galaxies, many numerical simulations 
have been performed and compared to observat ions (e.g. 
iBarnes fc Hernquisd Il99i iBurkert fc Naabl l200l and ref- 
erences therein). The general consensus is that it is pos- 
sible to generate remnants which agree in their kinematic 
properties with observed elliptical galaxies. Recent merger 
simul ations use self-consiste nt cosmological initial condition 
llKhochfar fc Burker t 2006), in simulations with only stars 
fe.g. lNaab fc Burker«200oTl and in simulation s with gas, star 
formation and feedback (e.g. lCox et aljEo05h . Usually these 
simulations neglect the past history of a possible progenitor 
galaxy, with the simulations starting from idealised initial 
conditions. 

Observationally, there is evidence that the initial con- 
dition s of merging galaxies cover a wide r ange, from being 
'dry' Jvan Dokkuml 120051: IBell et al.ll2005l) to being ultra- 
luminous (for a review see lsThweizei 127)051) . It is therefore 
expected that the stars ending up in a spheroid could not all 
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have been formed in the same kind of event, but must have 
formed in different events which depend on the past history 
of the merging galaxies. Since merger simulations indicate 
that violent relaxation is not complete and that the stars 
manage to reme mber their in itial energy and orbital angu- 
lar momentum l|Barnesl 119981. it is necessary to take into 
account the past history of a galaxy. The main task now is 
to connect these past histories with the properties of merger 
remnants. 

Large observational studies that collect data from sev- 
eral thousands of galaxies, like the Sloan Digital Sky Sur- 
vey, reveal that the galaxy population follows remarkable 
trends: e.g. the surface mass density of galaxies increases 
with mass until a characteristic ma ss scale of Mq = 3 x 10 1 
Mq at which it becomes constant (iKauffmann et aH f2003). 
The constant surface mass density is mainly associated with 
elliptical galaxies and galaxies having significant bulges. Fur- 
thermore, detailed studies of the size distribution of elliptical 
galaxies reveal that the scatter in sizes of elliptical galax- 
ies of a given mass is log-normal distribu ted with a scatte r 
which decreases for larger galaxy masses (|Shen et alj [20031. 
Those authors could explain the evolution of the sizes by as- 
suming continued mergers of galaxies which had initially all 
the same mass but sizes that followed an appropriate distri- 
bution. The origin of this distribution however remains un- 
solved. Another interesting observed correlation is that the 
surface mass density of g alaxies in different envi ronments 
only differs at low masses l|Kauffmann et al"ll2004Fl . 

In this paper, we address the question of where stars 
that end up in spheroids were formed. The paper is struc- 
tured as follows: we begin by explaining the model ingre- 
dients we use, followed by a section on the past merging 
history of galaxies. Then we introduce our definition of the 
two stellar components found in spheroids and discuss how 
these components evolve in the simulated galaxy population. 
In section five we present our conclusions. 



2 GALAXY FORMATION MODELLING 

The main strategy behind the modelling approach we follow 
is first to calculate the collapse and merging history of indi- 
vidual dark matter halos, which is governed purely by grav- 
itational interactions, and secondly to calculate the more 
complex physics of the baryons inside these dark matter ha- 
los, including e.g. radiative cooling of the gas, star formation, 
and feedback from supernovae by simplified prescriptions on 
top of the dark matter evolution. Each of the dark mat- 
ter halos will consist of three main components which are 
distributed among individual galaxies inside them: a stel- 
lar, cold, and hot gas component, where the latter is only 
attributed to central galaxies, which are the most massive 
galaxies inside individual halos. In the following sections, 
we will describe briefly the recipes used to calculate these 
different compo nents which are mainly based on recipes 
presen ted in e.g. ICole et alJ dl994 l2000l) IKauffmann et alJ 
i|l999l) (hereafter. K99) and lSpringel et alJ fcOOll) (hereafter. 
SOI), and we refer readers for more details on model imple- 
mentations to their work and references therein. 

Throughout this paper we use the following set of cos- 
mological parameters: Qo — 0.3, Qa = 0.7, fi(,/fio = 0.15, 
as = 0.9 and h = 0.65. 
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Figure 1. Average mass fraction of disc stars populating 
spheroids depending on the mass resolution M m i„ of the merger 
tree. 

2.1 Dark Matter Evolution 

We calculate the merging history of d ark matter halos ac- 
cording to the prescription presented in lSomerville fc Kolattl 
Jl999ah . This approach has been shown to produce merg- 
ing histories and progenitor distributions in reasonable 
agreement with results from N-body simulations of cold 
dark matter structure formation in a cosmological context 
JSomerville et all2000l) . The merging history of dark matter 
halos is reconstructed by breaking each halo up into pro- 
genitors above a limiting minimum progenitor mass M m i„. 
This mass cut needs to be chosen carefully as it ensures that 
the right galaxy population and merging histories are pro- 
duced within the model. Progenitor halos with masses below 
M m i n are declared as accretion events and their histories 
are not followed further back in time. Progenitors labelled 
as accretion events should ideally not host any significant 
galaxies in them and be composed mainly of primordial hot 
gas at the progenitor halo's virial temperature. The mass 
scale at which this is the case can in principle be estimated 
from the prescriptions of supernova feedback and reioniza- 
tion presented in section ?2. 2. II However, to achieve a good 
compromise between accuracy and computational time, we 
instead estimated M m i„ by running several simulations with 
different resolutions and chose the resolution for which re- 
sults in the galaxy mass range of interest are independent 
of the specific choice of M m i n . Changing the mass resolu- 
tion mainly affects our results at low galaxy mass scales as 
shown in Fig.l, leaving massive galaxies nearly unaffected. 
Throughout this paper we will use M mln = 2 x 10 9 Mq 
which produces numerically stable results for galaxies with 
stellar masses M, > 10 9 Mq . 

2.2 Baryonic Physics 

As mentioned above, once the merging history of the dark 
matter component has been calculated, it is possible to fol- 
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low the evolution of the baryonic content in these halos for- 
ward in time. We assume each halo consists of three com- 
ponents: hot gas, cold gas and stars, where the latter two 
components can be distributed among individual galaxies 
inside a single dark matter halo. The stellar components of 
each galaxy are additional divided into bulge and disc, to 
allow morphological classifications of model galaxies. In the 
following, we describe how the evolution of each component 
is calculated. 

2.2.1 Gas Cooling & Reionization 

Each branch of the merger tree starts at a progenitor mass 
of M m in and ends at a redshift of z = 0. Initially, each halo 
is occupied by hot primordial gas which was captured in the 
potential well of the halo and shock heated to its virial tem- 
perature T V ir = 35.9 [14/ ( km s" 1 )] K, wh ere V c is the cir- 
cular velocity of the halo (IWhite fc FrenkHl99ll. K99). Sub- 
sequently this hot gas component is allowed to radiatively 
cool and settles down into a rotationally supported gas disc 
at the cent re of the halo, which we identify as the central 
galax y (e.g Isilk1ll977t IWhite fc Reedll978t1white fc Freiikl 
1991). The rate at which hot gas cools down is estimated 
by calculating the cooling rad ius inside the halo us i ng the 
cooling functions provided bv ISutherland fe Dopital (Il993l) 
and the prescription in SOI. In the case of a merger between 
halos we assume that all of the hot gas present in the pro- 
genitors gets shock heated to the virial temperature of the 
remnant halo, and that gas can only cool down onto the 
new central galaxy which is the central galaxy of the most 
massive progenitor halo. The central galaxy of the less mas- 
sive halo will become a satellite galaxy orbiting inside the 
remnant halo. In this way, a halo can host multiple satellite 
galaxies, depending on the merging history of the halo, but 
will always only host one central galaxy onto which gas can 
cool. The cold gas content in satellite galaxies is given by 
the amount present when they first became satellite galaxies 
and does not increase, instead it decreases due to ongoing 
star formation and supernova feedback. 

In the simplified picture adopted above, the amount 
of gas available to cool down is only li mited by the 
unive rsal baryon fraction Q,i,h 2 = 0.024 iSpergel et alJ 
2003). However, in the presence of a photoionising back- 
grou nd the fraction of baryons captured in halos is reduced 
(e.g. [Navarro fc SteinmetJl997l :|Gnedin|2000}|B enson et alJ 
2002) and we use the recipe of ISomervillel j2002l) . which is 
based on a fitt ing formulae d erived from hydrodynamical 
simulations by iGnedinl ((2000) , to estimate the amount of 
baryons in each halo. For the epoch of reionisation, we as- 
sume Zreian = 20, which is in agreement with observations 
of the temperature-polarisatio n corre l ation of the cosmic mi- 
crowave background bv lKoeut et alJ i2003l) . 

2.2.2 Star formation in Discs and Supernova Feedback 

Once cooled gas has settled down in a disc, we allow for 
fragmentation and subsequent star formation a ccording to 
a par ameterised global Schmidt-Kennicutt law iKennicutd 
1998) of the form M, = aM co i d /t dyntga i, where a is a free 
parameter describing the efficiency of the conversion of cold 
gas into stars, and t d yn, g ai is assumed to be the dynamical 



time of the galaxy and is approximated to be 0.1 times the 
dynamical time of the dark matter halo (K99). As in K99 
we allow star formation only in halos of V c < 350 km/s to 
avoid too bright central galaxies in clusters. 

Feedback from supernovae plays an important role in 
regulating star formation in small mass halos and in pre- 
venting too massive satellite galaxies from forming. We im- 
plement feedback based on the prescription presented in K99 
with 

AMreheat = tn^LAM,. (1) 

Here we introduce a second free parameter e which repre- 
sents our lack of knowledge on the efficiency with which 
the energy from supernovae is going to reheat the cold gas. 
The expected number of supernovae per solar mass of stars 
formed is given by »7sjv = 5 x 10~ 3 ; take n as the value for 
the Scalo initial mass function iScaloll986t) . and Esn = 10 51 
erg is the energy output from each supernova. We take V c 
as the circular velocity of the halo in which the galaxy was 
last present as a central galaxy. 

2.2.3 Galaxy Mergers 

We allow for mergers between galaxies residing in a single 
halo. As mentioned earlier, each halo is occupied by one 
central galaxy and a number of satellite galaxies depending 
on the past merging history of the halo. Whenever two halos 
merge, the galaxies inside them are going to merge on a time- 
scale which we calculate by estimating the time it would take 
the satellite to reach the centre of the halo under the effects 
of dynamical friction. Satellites are assumed to merge only 
with central galaxies and we set up their orbits in the halo 
according to the prescription of K99, modified to use the 
Coulomb logarithm approximation of S01. 

If the mass ratio between the two merging galaxies 
is M ga i tl /M gaL2 < 3.5 (Mgai,i ^ M ga i a ) we declare the 
event as a major merger and the remnant will be an ellip- 
tical galaxy and the stellar components and the gas will be 
treated according to the prescriptions below. In the case of 
minor merger M ga i,i/M ga i,2 > 3.5 the cold gas in the disc 
of the smaller progenitor is assumed to settle down in the 
gas disc of the remnant and its stars contribute to the bulge 
component of the remnant (e.g. K99). 

2.2.4 Formation of Ellipticals and Bulges 

iToomre fc Toomrd il972f) suggested that major mergers will 
lead to the formation of elliptical galaxies. Indeed detailed 
numerical simul ations in the last decade seem to support this 
hypo thesis (e.g. iBarnes fc Hernauisdfl992l : iBurkert fc Naabl 
2003, and reference therein), and we will assume in the fol- 
lowing that major mergers disrupt the discs in progenitor 
galaxies, as seen in various numerical simulations, and relax 
to a spheroidal distribution. During the merger, any cold 
gas in the discs of the progenitor galaxies is assumed to be 
funnelled into the centre of the remnant where it ignites 
a starburst which transforms all of the cold gas into stars 
contributing to the spheroidal component (e.g. K99; S01; 
and references therein). The second assumption is certainly 
a simplification of what might happen since we neglect the 
possibility that not all of the cold gas is funnelled to the 
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Figure 2. Predicted /-band Tully-Fisher relation for Sb/Sc 
galaxies in our simulation. The solid line is the fit to the observed 
data bv iGiovanelli et alj Jl997l) and dashed lines give their esti- 
mate for the scatter. 



Figure 3. The predicted i\-band luminosit y function in our 
sim ulation in compariso n to the observations bv lCole et alj J200l[) 
and lHuang et alj l20o3l . 



centre but some fraction of it can e.g. settle down into an 
extended disk which continues g rowing inside-out by fresh 
supply of gas from tidal tails (e.g. lBarne s fc Hernauistl 19911 
i Mihos fc Hernauistll99"5lBarnesll200ll.l2002l) . The results of 
iBarnesI feOOSfl indicate that 40% - 80% of the initial gas 
in the disc could end up in the central region of the rem- 
nant and be consumed in a starburst. The exact amount is 
somewhat dependent on the merger geometry and on the 
mass ratio of the merger. Unfortunately, a large survey in- 
vestigating the gas inflow to the centres of merger remnants 
is not available to date so that we use the simplified ap- 
proach of assuming that all cold gas gets used up in the 
central starburst. This prescription for the fate of the cold 
gas results in an overestimate of the spheroid masses and 
an underestimate of the secondary disc components in our 
model. This is not very significant for massive ellipticals 
since the y are mainly formed in rel atively weak dissipative 
mergers (iKhochfar fc Burkerdl2003l) . 

Another simplifying assumption is that we neglect the 
feeding of super massive black holes in the centre of the 
remnant or feedb ack effects on the gas f rom t he central 
source. However, lHaehnelt fc Kauffmannl (2000) estimate 
that a cold gas mass fraction of less than 1% accreted onto 
the black hole is sufficient to recover the M, — a relation 
and we therefore neglect this effect on the amount of gas 
available for the central starburst. A larger effect on our re- 
sults arises from the fact that we neglect feedback from the 
central source into the surrounding intergalactic medium. 
As a consequence, our estimates of stars formed in a central 
starburst will be too high and should be viewed as upper 
limits. 



2.2.5 Free parameters and predicted present-day 
properties of modelled galaxies 

We normalise our two model parameters for the star forma- 
tion efficiency a and feedback efficiency t by matchin g the 
/—band Tully-Fisher relation of lGiovanelli et aD (ll997tT and 
requiring that spiral central galaxies of halos with circular 
velocity Vc = 220 km/s have on average 10 11 Mr of stars 
and a few times 10 9 Mq of cold gas iSomerville fc Kolatd 
1999b). The best fit model parameters we get in this way 
are e = 0.2 and a — 0.1. 

Figure shows our predicted I-band Tully-Fisher rela- 
tion for Sb/Sc galaxies in comparison to the observed one by 
IGiovanel li et alJ (119971) . We selecte d Sb/Sc galaxies follow- 
ing lSimieri^^ieVaucouleursl (Il986l) according to their bulge 
to total light in the range 1.5 < M B ,bui ge — Mb, total < 2.2. 

The model predictions for the 7^-band luminosities of 
the overall galaxy population and the early-type galaxy lu- 
minosity are shown in Fig. 21 & refmod2, respectively. We 
choose early-type galaxies to have bulge to total light ra- 
tios with 1 ^ M b bulge — Mb. total, which includes SO galax- 
ies according to ISimien fc de Vaucouleursl (Il986l) . Our lu- 
minosity functions somehow over predicts the abundance 
of very luminous objec ts, which is a know n problem and 
previously discussed in iBenson et al.l <l2003l) . The solution 
to this proble m is commonly believed to come from AGN- 
feedback. ICroton et alJ l|200tj) who use semi-analytic mod- 
elling based on the same recipes that we do find a luminosity 
function very similar to ours (see Fig. 8 of their paper) . How- 
ever, the inclusion of AGN-feedback in their models reduces 
the amount of very luminous galaxies and fits the data very 
well at the high mass end. As for our results presented in 
this paper the neglect of AGN-feedback results in a few to 
bright and too massive early-type galaxies as can be seen in 
Fig. |21 However, these galaxies are rare and extreme cases 
not considered in our analyses. Our early-type luminosity 
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Figure 4. The predicted K-hand luminosity function in our sim- 
ul ation for early-type gala xies in comparison to the observations 
bv lKochanek et alj 1200 ll) . 



fiction fits the data very well over most of the considered 
range and we therefore delay further discussion of AGN- 
feedback to a forthcoming paper (Khochfar & Silk in prep.). 



3 PAST MERGING HISTORY OF 
PROGENITOR GALAXIES 

The composition of stars in bulges and elliptical galaxies 
is very dependent on the previous merging history of the 
progenitors and their ability to cool gas and regrow large 
stellar discs in between major mergers iKhochfar fc Burkerd 
( 2003) showed that based on the stellar mass of the rem- 
nant, a transition appears at which regrowing of a stellar 
disc becomes inefficient. This leads to major mergers mainly 
between bulge-dominated galaxies with too low fractions of 
available cold gas to ignite central starbursts, sometimes 
called "dry" mergers. It appears that such a transition is 
necessary to explain propert ies like e.g. the isophotal sh ape 
of massive elliptical galaxies jKhochfar fc Burkerdfeooit) . To 
stress the importance of taking into account the past merg- 
ing history of progenitors, we investigate in Fig.[5]the ques- 
tion of the relative number of mergers between galaxies with- 
out bulges to the number of mergers between galaxies with 
bulges, depending on the masses of the remnants. As can be 
seen, it becomes very unlikely to find mergers between two 
pure disc galaxies that lead to very massive remnants. At 
masses above M* ~ 3 x 10 9 M©, the majority of mergers 
take place between galaxies that have already experienced 
a major merger in their past. It is interesting to note that 
we do not find any mergers between bulge-less disc galaxies 
in our simulation at masses larger than the characteristic 
mass scale Mc, suggesting that the effect of constant sur- 
face mass density /x* in elliptical galaxies more massive than 
Mc is closely related to the bulge fractions in the progen- 



Figure 5. Ratio of the number of mergers between galaxies with 
bulges and without bulges as a function of the remnant mass 
found in our simulations. Th e solid line shows the results for the 
cosmological mass function o fl Jenkins et alJfcOOll) . and the other 
lines examine results for different environments. The vertical dot- 
dashed line ma rks the characteristic ma ss scale Mc = 3 X 10 10 
Mq reported in lKauffmann et alj l2003t) 



itor galaxies. Inspecting the environmental dependence, we 
find that the ratios evolve independently of the environment 
of the galaxy and that bulge-less mergers stop occurring at 
lower masses in dense environments like groups and clus- 
ters. This is a consequence of the hig her merger fractions 
at earlier times in these environments jKhochfar fc Burkerd 
1200 ID. 



4 MERGER AND QUIESCENT 
COMPONENTS IN SPHEROIDS 

The populations of stars present in a spheroid at a given 
time in its evolution are a composite of stars formed in var- 
ious progenitors at different times and in different modes. 
Here we distinguish between two modes, the merger-induced 
'merger component' and the disc mode 'quiescent compo- 
nent'. The merger component is a composite of all stars that 
were formed because of the consumption of cold gas during 
major mergers in progenitor galaxies. On the other hand, 
the quiescent component is the amount of stars formed in 
gaseous discs according to the Schmidt-Kennicutt law dur- 
ing the evolution of the progenitor galaxies. When calcu- 
lating each component for a spheroid, we sum up the mass 
of the merger component so far and subtract it from the 
spheroid mass to get the quiescent component. 

To investigate systematic trends in the fraction of 
merger-induced and quiescent stellar components in ellip- 
tical galaxies and bulges, we use in the following the condi- 
tional distribution 

4>(M q /M bul ,M*) 



p(M q /M bu i\M*) 



(2) 



^{M'JM' buV M,)d{M' q /M' bui y 
where M q is the stellar mass in the spheroid previously 
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Figure 6. Conditional distribution of the fraction of quiescent 
stars ending up in spheroids as a function of the total galaxy 
mass for the overall galaxy population at redshift z $C 0.3. The 
solid line shows the mean of the distribution and the vertical dot- 
dashed line marks the characteristic mass scale Mc- 

formed quiescently, Mbul the mass of the spheroid and M* 
the total stellar mass of the galaxy. The merger component 
in the spheroid is then simply given by 1 — Mq/Mbul- If not 
stated otherwise, this distribution is derived for the simu- 
lated galaxy population at redshifts z ^ 0.3. 

The quiescent component of stars in spheroids is get- 
ting more dominant with galaxy mass up to masses around 
Mc when it starts becoming a constant fraction of the stars 
in the spheroid (Fig. The behaviour at galaxy masses 
below Mc can be understood by acknowledging that mas- 
sive spheroids form late in the hierarchical galaxy formation 
paradigm, thereby allowing more stars to be formed in pro- 
genitor discs between individual major merger events, and 
that feedback from supernovae is less efficient in reheating 
cold gas in discs which are embedded in massive halos. It 
is worth noting that even spheroids in galaxies as small as 
M* = 10 9 Mq have on average only ~ 30% of their stars 
being formed in past merger-triggered central starbursts, in- 
dicating that the main mode of star formation is taking place 
in discs. The quiescent fraction of stars can get as large as 
~ 85% at masses above Mc, where it then becomes con- 
stant. 

4.1 Redshift Evolution 

The distribution found for the galaxy population at z ^ 
0.3 is likely to be dependent on redshift. Galaxies at high 
redshifts tend to have a higher fraction of cold gas available 
in their discs than their counterparts at low redshift when 
they participate in mergers as shown in Fig.|5] This is mainly 
due to the time available to make stars in the disc before the 
major merger and the very short cooling times for the hot 
gas. To investigate the build up of the distribution at low 
redshift, we look at the trajectory of a galaxy defined by its 
positions in the M q /Mbui — M* plane through time. Since 



Figure 7. Trajectory of a representative galaxy with final mass 
M* = 1.5 X 10 11 Mq at z = 0. We show the evolution as a function 
of the bulge mass Mbul and as a function of the total stellar mass 
M„, indicated by the symbols and the solid line, respectively. The 
average of the overall galaxy distribution is displayed by the dot- 
ted line, and the vertical dot-dashed line marks the characteristic 
mass scale Mc- Please note that the trajectory as a function of 
total stellar mass starts around M, ~ 10 9 Mq because we start 
following trajectories after the first occurrence of a major merger 
between progenitor galaxies. For guidance wc label the bulge and 
total stellar mass at the first occurrence of a significant central 
star burst with the corresponding redshift of that event (z = 4.42) 
and at a redshift of z = 0. 

the present day value of M q is integrated over the entire 
history of a galaxy, we calculate the trajectories by summing 
up at each redshift the quiescent stellar components in the 
individual progenitor bulges and ellipticals and divide it by 
the sum of the stellar mass of the progenitor spheroids: 

Mbul(z) Y,i=x,N prog M ^l,i{ z Y 

Note that elliptical galaxies might contain a small disc which 
is in the process of regrowing, and that we do not attribute 
these stars to the spheroid when calculating the trajectory. 
A typical trajectory for an elliptical galaxy more massive 
than Mc is presented in Fig. Q This galaxy experienced a 
gas-rich major merger at a redshift of z = 4.42 between two 
progenitors leading to a significant decrease in the quiescent 
fraction Mq/Mbul- The significant decrease is not only due 
to the amount of cold gas available in the merger but also to 
the fact that most of the stars in the individual progenitors 
are still in their discs, and that only a few, if any, galaxies 
consist of spheroids, which makes a single major merger have 
a large impact. Once significant spheroids in the population 
of progenitors start to exist, major mergers between two 
progenitors tend to lower the overall fraction of quiescent 
stellar components only slightly. 

The main mode by which M q /Mt u i increases is by mi- 
nor merging with satellite galaxies. These satellite galax- 
ies, though smaller in mass than their central galaxies, have 
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Figure 8. The evolution of the mean value Mq/M^ul of the con- 
ditional distribution p(M q /Mf, u i\M t ) as a function of stellar mass 
M, in five different redshift bins between 2 = and 2 = 4. The 
vertical dot-dashed line marks the characteristic mass scale Mq . 



Figure 9. The evolution of the mean gas fraction in progenitors 
of major mergers as a function of stellar mass M* in five different 
redshift bins between 2 = and 2 = 4. Please note that M* is 
the stellar mass of the merger remnant. 



grown stellar discs which are massive enough to increase 
M q /Mi, u i. That minor mergers are important becomes evi- 
dent by noting that they are around an order of magnitude 
more frequent than major mergers. Once the bulge mass 
reaches roughly Mc the ratio M q /Mt u i is around the av- 
erage value inferred from the overall distribution at that 
redshift. Subsequently the mergers, minor as well as major, 
only create a small scatter around this value, consistent with 
the assumption of bulge-dominated mergers. 

The trajectory in Fig. Q suggests that the mean value 
of the conditional distribution will change with time, in 
co-evolution with the merger activity. Fig. |H] compares the 
mean value of Mq/M^ui found in the conditional distribu- 
tions p(M q /Mbui\M t ) in five different redshift bins. At ear- 
lier times, the means show systematic lower fractions of qui- 
escent components in their spheroids. This is connected to 
the higher fraction of available gas in major mergers at ear- 
lier times as can be seen in Fig. |5] The different averages 
show similar behaviours with stellar mass below Mc only 
differing by constant offsets and a slightly steeper increase 
of M q /Mi, u i toward Mc- Above Mc, the starburst compo- 
nent in spheroids starts to increase again by up to ~ 30% 
depending on the redshift, a result of massive galaxies at 
high redshift likely being formed in a very gas-rich merger 
event, or having bulges that were formed very early on. At 
this point, it is not clear if this trend at high redshifts is 
real or just a product of not including feedback from super- 
massive black holes during major mergers and we will defer 
discussion of possible effects of AGN-feedback to a following 
paper (Khochfar & Silk in prep.). However, note that we 
usually do not find more than a few galaxies in the highest 
mass bins, which makes finding such galaxies very unlikely. 
From the different curves in Fig. [SI it is not obvious whether 
Mc plays the same fundamental role at all redshifts or if the 
characteristic mass scale changes with redshift. In any case, 



the differences according to our simulations are not very 
large for the overall galaxy population. 

To understand the shape of the curves in Fig. |H| it is 
important to understand the extent to which the evolution 
in Mq/Mtui at each redshift and mass scale is driven by 
galaxies with different morphologies. In Fig. HUI we compare 
the overall stellar mass in spheroids divided by the total 
stellar mass present at different redshifts and masses. The 
closer to Mbui/M* = 1, the more likely it is that most of the 
galaxies at a given mass scale are elliptical galaxies which 
have recently formed and did not regrow significantly large 
discs since then. At each redshift, our model predicts that 
the most massive galaxies, which are less massive with in- 
creasing redshift, are recently formed elliptical galaxies. As 
a consequence, the quiescent component in these galaxies 
resembles the average composition that results from major 
mergers at this redshift. Going to smaller galaxies, the bulge 
mass fractions decrease, indicating that more mass is found 
in stellar discs, which in turn means that the bulges must 
have formed earlier to allow for the growth of large stel- 
lar discs. As argued above, early mergers are more gas-rich 
and therefore lead to smaller ratios of M q /Alb u i- This, and 
the fact that galaxies with increasing mass have a larger 
number of minor mergers, which in general increase the qui- 
escent component of bulges, accounts for the increase in the 
quiescent component in spheroids with increasing mass. 

It is interesting to note that Fig. 1101 suggests that 
the characteristic mass scale Mc approximately marks the 
transition point at which most of the stars start to be in 
spheroids independent of the redshift. Below Mc, the av- 
erage distribution of stars in bulges and discs is roughly 
independent of the redshift, and shows the same mass de- 
pendence. Only at galaxy masses above Mc does the distri- 
bution of stars in spheroids and discs become significantly 
redshift-dependent. This behaviour can be attributed to two 
features of the CDM paradigm: first, the decline in the 
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Figure 10. Bulge to total stellar mass fraction in five different 
redshift bins between 2 = and 2 = 4. We calculated the fraction 
by summing up all stars present in spheroids at a given redshift 
and mass scale and dividing this value by the total amount of 
stars available at this mass scale. The vertical dot-dashed line 
marks the characteristic mass scale Mc- 



Figure 11. The evolution of the mean value Mq/M^ul of the con- 
ditional distribution p(M q /Mt, u i \M*) as a function of stellar mass 
M„ in three different environments. We defined environments by 
the the mass of the dark halo the galaxies reside in. Halos with 
10 12 M Q , 10 13 M Q and 10 14 Mq are classified as field, small 
group and group/massive group environments, respectively. 



merger rate going to lower redshifts and the connected abil- 
ity of stellar discs to grow longer without being disrupted, 
and second, most massive galaxies at a given redshift are 
likely to be those which formed in a major merger, resulting 
in an elliptical galaxy. 

To summarise, the build up of bulges and elliptical 
galaxies in a merger-driven model suggests that early gas- 
rich mergers lead to progenitors with a high fraction of stars 
in their spheroids, and which were formed during a dissi- 
pative central starburst event. This fraction is significantly 
higher than that found at later times in spheroids. While the 
progenitors of a present-day galaxy continue to grow stellar 
discs, minor mergers occur. These events tend to increase 
the quiescent fraction of stars in spheroids by contribut- 
ing their disc stars to the spheroid of the remnant. These 
minor mergers occur mainly during the peak times in the 
merger activity which is dependent on the environ ment that 
the galaxy resides in (iKhochfar fc BurkerdEool . As soon 
as the bulge mass of the progenitor galaxies reaches ~ Mc, 
only small changes occur. However, depending on the mor- 
phology of the resultant galaxy the final stages might differ 
as well as the characteristic mass scale (Khochfar &Silk, in 
preparation). 

4.2 Environmental Dependence 

In paragraph s. II we presented the evolution of each individ- 
ual stellar component in spheroids. This evolution is very 
dependent on the merger activity that galaxies experienced 
in their past. It is to be expected that the different evo- 
lutions of the merger rates with redshift in field and high 
density environments leave a characteristic signature in the 
merger component of spheroids. 

In Fig. 1111 we compare the fraction of quiescent com- 



ponents in the spheroids of galaxies found in different envi- 
ronments in our simulation. We classify environments into 
three different categories according to the mass of the dark 
matter halo the galaxies reside in. Galaxies in dark matter 
halos of 10 12 M Q , 10 13 M and 10 14 M are classified as 
field, small group, and group/massive group environments 
in our simulation. The average fraction of quiescent com- 
ponents in the bulges of spheroids less massive than ~ Mc 
is systematically higher in field environments meaning their 
merger components are smaller. Furthermore we find in our 
simulations that there is no systematic difference for galaxies 
above ~ Mc- The difference between environments dimin- 
ishes, the closer the galaxy mass is to Mc- Recalling the 
results of section it becomes clear that massive galaxies 
in different environments are unaffected because the contin- 
uous growing by mergers, minor as well as major, erases the 
memory of the first merger epoch in different environments. 
Only those galaxies significantly less massive than ~ Mc did 
not grow by many mergers after their initial mergers, but 
survived until today without much interactions as satellite 
galaxies. Since the initial merger is very much environmental 
dependent, gas-rich early mergers in high density environ- 
ments and less gas-rich late mergers in field environments, 
differences arises between environments. 



5 DISCUSSION AND CONCLUSION 

In this paper we have investigated the origin of stars ending 
up in bulges and elliptical galaxies for systematic trends with 
mass, redshift and time of the last major merger. We divide 
stars in spheroids into two categories, merger and quiescent, 
according to their formation mode. 

The merger component is associated with major merg- 
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ers in which cold gas in progenitor discs is driven to the 
centre of the remnant and a central starburst is ignited. 
Our simulations predict that the fraction of the merger com- 
ponent in spheroids is an increasing function of formation 
redshift and a decreasing function of spheroid mass. The 
latter is due to the longer time that progenitors of massive 
galaxies have to grow discs, and an enhanced number of 
minor mergers, which generally decreases the fraction. The 
higher fraction at earlier formation times is a result of the 
short cooling time scales for gas and the insufficient build 
up of stellar discs before the final major merger. Our sim- 
ulations predict that the average merger component in the 
local galaxy population can be as high as 30% of the bulge 
mass in galaxies with M* ~ 10 9 Mq and as low as 15% in 
galaxies with M, > M c = 3 x 10 10 M©. 

The quiescent component is associated with stars pre- 
viously formed in gaseous discs according to the Schmidt- 
Kennicutt law and later being added to the spheroid during 
major mergers. We find that the fraction of quiescent stars 
dominates bulges and elliptical galaxies. For galaxies above 
the critical mass scale, the fraction of quiescent stars be- 
comes constant at a value of M q /M^ u i ~ 85%, the result 
of mergers, minor as well as major, only introducing small 
deviations as soon as galaxies are more massive than Mc- 

The fraction of stars in bulges and discs shows a re- 
markable relation with the galaxy mass scale. At galaxy 
masses below Mc , we find only a weak redshift dependency 
and a decreasing fraction of stars in bulges with decreasing 
mass. However, at masses above Mc we find a strong red- 
shift dependence. At higher redshifts, the fraction of stars 
in bulges becomes larger in massive galaxies. This means 
that the most massive galaxies around at each redshift are 
likely to be elliptical galaxies and that only at later times 
do massive spiral galaxies start occurring. 

Furthermore, our simulations indicate that mergers 
leading to remnants with mass larger than Mc always in- 
clude progenitors with bulges and that more than 50% of 
the mass in progenitors of mass M* ~ Mc is in bulges. This 
has an important influence on the role of AGN-feedback as 
these bulges will harbour super-massive black holes. 

Looking at the environmental dependence, we find that 
only galaxies less massive than Mc show differences. This is 
mainly related to the different epoch at which these galaxies 
assembled. Assembly is faster in high density environments 
resulting in larger merger components in bulges. Galaxies 
above Mc do not show any significant environmental de- 
pendency because of continued merger activity in contrast 
to smaller galaxies. 

Recently ISpringel fc Hernauistl i2005l) suggested that 
mergers between extremely gas-rich discs could lead to the 
formation of a disc galaxy during a major merger. The au- 
thors argued that these kind of mergers are likely to occur 
at high redshifts. However, our simulations indicate that the 
maximum gas fraction is at ~ 40% in mergers at high red- 
shift and most likely between 20% — 30% for a wide range of 
masses and redshifts. In general the gas fraction in mergers 
is a decreasing function of mass. At very large masses and at 
high redshifts a small population of gas rich massive merg- 
ers occurs, which tends to slightly deviate from this trend. A 
reason for these events in our simulation could be associated 
with missing feedback from AGNs. Those massive gas-rich 



mergers at high redshift are likely to have progenitors with 
black holes. 

Simulations bv lSpringel fc Hernauistl i2005T) show that 
in general the merger component is more centrally concen- 
trated than the quiescent component. The effective radius 
of the merger component in their simulation is ~ 5.7 times 
smaller than that of the quiescent component. This suggest 
an interesting behaviour for the sizes of two galaxies of the 
same mass but different merger component. The remnant 
from the merger including less gas and hence less merger 
component will have a larger sizes. As we have shown above 
remnants at higher redshifts have larger merger fractions. 
This and the observational indicatio n that disc galaxies 
at higher redshifts are smaller (e.g iTruiillo et all 120051) 
might be able to account for the ob s erved size-evolution of 
elliptical galaxies JPaddi et al.ll2005l : ITruiillo et alJl2005f) . 

We would like to thank Ignacio Trujillo and Emanuele 
Daddi for pointing out the high redshift data. SK acknowl- 
edges funding by the PPARC Theoretical Cosmology Rolling 
Grant. 
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